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We discuss the main features and predictions of the GLMM model, which is based on a 
QCD motivated theoretical approach, and successfully describes the experimental data on 
total, elastic and diffractive cross sections. In addition we calculate the survival probability 
for a SM Higgs at the LHC, and compare our results with those of the Durham group. 



1 Introduction 

Over the past few years the subject of "soft physics" has reemerged from the shadows, and 
has aroused the interest of the phenomcnological community. This in no small way, due to 
the realization that the calculation of the probability of detecting a diffractive hard pQCD 
process e.g. Higgs production at the LHC, also depends on the underlying secondaries which 
are produced by "soft" rescattering. Central diffractive production e.g. (Higgs boson, 2 jets, 

2 7's, Xc), are accompanied by gaps in rapidity, between the two outgoing projectiles, and the 
centrally produced particles, which makes their detection easier. The subject of the survival 
of these rapidity gaps was initiated over twenty years ago [T], and has been refined over the 
interim period [2] , [5] . 



2 Details of our Model 

The details on which our two channel (GLM) model is based i.e. the Good- Walker (G-W) 
mechanism U can be found in [2J. See also U. Maor's talk in these proceeedings. It is well 
known that G-W neglects the diffractive production of large mass states (Mueller diagrams [5]), 
and to successfully describe the diffractive data, these need to be included by adding the relevant 
triple Pomeron contributions (see L.H. diagrams in Fig.l). In addition it is also necessary to 
include diagrams containing Pomeron loops, as shown in the R.H. diagrams in Fig.l. 

To simplify the problem of summing the Pomeron loop diagrams, we assume that at high 
energies only the triple Pomeron interaction is essential, this conjecture has been proved in 
peturbative QCD (see Refs.[[6] [7]]). Mueller [8] has shown that in the leading log x approx- 
imation of pQCD for a large number of colours (N c 3> 1), the correct degrees of freedom arc 
colourless dipoles. 
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2.1 Summing Interacting Pomeron Diagrams 

In the leading order approximation of pQCD, only one Pomeron (dipole) splitting into two 
Pomerons (dipolcs), and two Pomerons (dipolcs) merging into one Pomeron (dipole) are con- 
sidered. All other Pomeron vertices do not appear in the leading log x approximation of pQCD. 
We therefore restrict ourselves to sum only Pomeron diagrams containing triple Pomeron ver- 
tices. We add a caveat, that we neglect the 4JP term (which is needed for s channel unitarity), 
however, this term is only significant at energies W > 10 5 , which is at the limit of the validity 
of our model. 

To make the calculation tractable we further assume that the slope of the Pomeron trajectory 
a'p = 0. The results of our numerical fit to the relevant data, which we will discuss later, 
(a'p = 0.01) lends credence to this assumption. 

The theory which includes all the above ingredients can be formulated in terms of a gener- 
ating function [9l [TO] 

Z(y, u) = ]T P„(y) u n , (1) 

n 

where, P n (y) is the probability to find n-Pomerons (dipoles) at rapidity y. At rapidity y = 
Y = ln(s/so) we can impose an arbitrary initial condition. For example, demanding that there 
is only one fastest parton (dipole), which is Pi(y = Y) = 1, while P„>i(y = Y) = 0. In this 
case we have the following initial condition for the generating function 

Z{y = Y)=u. (2) 

At u = 1 

Z(y, u = 1) = 1, (3) 

which follows from the physical meaning of P n as a probability. The solution, with these two 
conditions, will give us the sum of enhanced diagrams. 

For the function Z (u) the following simple equation can be written (see Ref . [2] and references 
therein) 

- 9 -^ = -T{l^2)u{l-u) d -^l + r(2-l)u(l-u)^4, (4) 
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where, T(l — > 2) describes the decay of one Pomeron (dipole) into two Pomerons (dipoles), 
while r(2 — > 1) relates to the merging of two Pomerons (dipoles) into one Pomeron (dipole). 

Using the functional Z, we calculate the scattering amplitude [TUl [TT] . using the following 
formula: 

N(Y) ee ImA el (Y) = f^^T- ^ f^' ~ U=i ^(Y = Y , b), (5) 

; ri! a" u 

n— 1 

where, 7 n (F = io,b) is the scattering amplitude of n-partons (dipoles) at low energy. 

The generating function approach has the advantage that it can be solved analytically (see 
Ref. [12j). using the MPSI [T3] approximation. The exact expression for the Pomeron Green's 
function is given by 

G*<Y) = 1 " exp (-^) P (0,^) , (6) 

where r(0,a;) is the incomplete gamma function, and T (Y) = je AlpY . 7 denotes the am- 
plitude of the two dipole interaction at low energy. The MPSI approximation only takes into 
account the first term of the expression of the enhanced diagrams, neglecting other terms, as 
they are suppressed as e~ AY . Consequently, this approximation is only reliable in the region 
Y < mini-, -r±s]. 

3 Determining the Parameters of the Model and Results 
of the Fit 

The pertinent details of our fit to the experimental data, and our determination of the relevant 
parameters of the model, needed to describe the soft interactions, are contained in [2]. In 
this section we only mention the salient features, and results of the fit. Our fit is based on 
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Table 1: Fitted parameters for GLMM GW and GW+.P-enhanced models. 



55 experimental data points, which includes the p-p and p-p total cross sections, integrated 
clastic cross sections, integrated single and double diffraction cross sections, and the forward 
slope of the elastic cross section in the ISR-Tevatron energy range. The model gives a good 
reproduction of the data, with a \ 2 /d.o.f. m 1. In addition to the quantities contained in the 
data base, we obtain a good description of the CDF [14] differential elastic cross sections and 
the single diffractive mass distribution at t = 0.05 GeV 2 . An important advantage of our 
approach, is that the model provides a very good reproduction of the double diffractive (DD) 
data points. Other attempts to describe the DD data e.g. (see Refs.|15|.[3j). were not successful 
in reproducing the DD experimental results over the whole energy range. 

In Table 1 we list the values of the parameters obtained by a least squares fit to the experi- 
mental data, both for the G-W formalism (elastic data only), and for the G-W formalism plus 
enhanced graphs (elastic plus diffractive data). In Table 2 we compare our results with results 
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Figure 2: Energy dependence of GLMM cross sections. 



of the two versions of the model proposed by the Durham group [T5] and [3]. In Fig. 2 we display 
the predictions of the GLMM model's values for the various cross sections. Note that <7 e z and 
<j s d have completely different energy dependence, unlike the predictions of [3]. At an energy of 
7 Tev the predictions of GLMM are (in mb): (Jtot= 86.0, c e ;= 19.5, a s d= 10.7, add = 5.9 and 
the forward slope B ei = 19.4 GeV~ 2 . 





Tevatron 
GLMM KMR(07) KMR(08) 
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GLMM KMR(07) KMR(08) 


W=10 b GeV 
GLMM KMR(07) KMR(08) 
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11.8 13.3 19.0 


14.4 15.7 24.2 
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5.4 7.2 


6.1 13.4 
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0.43 0.46 


0.42 0.53 


0.41 0.57 



Table 2: Comparison of GLMM, KMR(07) and KMR(08) cross sections in mb. 




4 Survival Probability for Central Diffractive Production 
of the Higgs Boson 

A general review of survival probability calculations can be found in [15]. We denote by (| S 2 |) 
the probability that the Large Rapidity Gap (LRG) survives, and is not filled by secondaries 
from eikonal and enhanced rescattering effects (see Fig. 3). The expression for the survival 
probability can be written (| S\ ch |) = where 



N(s) 
and 



Jd 2 hd 2 b 2 E a < P \i> 2 <p\k> 2 A* H { s M)A k H { s M){i-Af{{ s ,{h l+ h 2 ))) 



D(s) = Jd 2 h d 2 b 2 [J2 l k < p\i > 2 < p\k > 2 A* H (s, h)A k H (s, b 2 ) 

As(b,s) denotes the "soft" strong amplitude of our model [2]. While for the "hard" am- 
plitude An{b, s), we assume an input Gaussian b dependence, i.e. Af k = Aff(s)rf k (b) 
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a) b) 

Figure 3: a) the survival probability in the G-W mechanism, b) illustrates the origin of the 
additional factor (| S 2 nh |). 

_ 2b 2 

with Tf k (b) = y e (Ri ' k> ■ The "hard" radii are constants determined from HERA 

data on elastic and inelastic J/vP production. We introduce two hard b-profilcs A p ^(b) = 

|g| exp (-5%) , and A p f f (b) = exp (-^) ■ The values 5^=5.0 GeV~ 2 and 

Bg=l GeV- 2 have been taken from ZEUS data. The value B%= (3.6) GeV~ 2 was used in [2J, 
this has now been changed in light of the latest measurements of the "hard" slope, by the HI 
group. This is in contrast to KMR treatment Q~5] where they assume: A p ^(b) = A P H '(b) oc 

exp 2S 73 ") w ith = Bf ne i = 4 or 5.5 GeV~ 2 The sensitivity of our results to the parame- 
ters of the "hard" amplitude are shown in Fig. 4 (the L.H. figure), note that for Bg =1 GeV~ 2 , 
changing the value of B^ from 3.6 to 5.0 GeV~ 2 , increases (| S 2 ch \) by m 70 % . Our results 
for (| S 2 |) = (| S\ ch |) x (| S 2 nh |) is given by the full line in Fig.4 (R.H. figure), it decreases 
with increasing energy, due to the behaviour of (| S 2 nh |). 

Our results and the Durham group's results for Survival Probability are given in Table 3. 
At an energy of 7 TeV we predict a value (| S 2 |) « 0.6%. We have also succeeded in summing 





Tevatron 
GLMM KMR(07) KMR(08) 


LHC (14 TeV) 
GLMM KMR(07) KMR(08) 


W=10 b GeV 
GLMM KMR(07) KMR(08) 


s'L h (%) 


v 5.3 2.7-4.8 


3.9 1.2-3.2 


3.2 0.9-2.5 




28.5 100 


6.3 100 33.3 


3.3 100 


S 2 {%) 


1.51 2.7-4.8 


0.24 1.2-3.2 1.5 


0.11 0.9-2.5 



Table 3: Comparison of results obtained for Survival Probability in Tel Aviv and Durham 
models 



the semi-enhanced contribution (see R.H. side of Fig.l) to the Survival Probability, and find 
that it is almost energy independent, and has a value ~ 100% at Tevatron and LHC energies 

5 Discussion and Conclusions 

We present a model for soft interactions having two components: (i) G-W mechanism for 
clastic and low mass diffractive scattering, and (ii) Pomeron enhanced contributions for high 
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Figure 4: L.H. figure: The dependence of S 2 at the LHC on and Bf n . R.H. figure: Energy 
dependence of S 2 for centrally produced Higgs. The full (dashed) line is for B^ = 5.0 (3.6) 
GeV~ 2 . 



mass diffractive production. In addition we find from our fit that the slope of the Pomeron 
a'p w 0.01. This is consistent with what one expects in pQCD, since for a BFKL Pomeron 
a'jp oc 1/Q 2 — > as s — > oo. Having a P — ► 0, provides a necessary condition that links 
strong (soft) interactions with the hard interactions described by pQCD. A key hypothesis in 
our model is that the soft processes are not "soft", but originate from short distances. We have 
only one Pomeron. There is no requirement for a "soft" and "hard" Pomeron. This is in accord 
with Hera data for Fa, which is smooth throughout the transition region [18] . 

To illustrate our achievements and problems, we compare our approach with the work of 
the Durham group |15j [3] . The main difference in the underlying philosophy of the two groups 
is that, the Durham approach is based on the parton model where there is only a short range 
rapidity interaction between partons, while we, due to exchange of gluons in QCD, have a long 
range rapidity interaction. Both approaches consider a'p as being small. In both programs 
the Pomeron interaction was taken into account. The difference between the two approaches is 
that KMR made an ad hoc "reasonable" assumption, that the multi-Pomeron vertices have the 
following form, for the transition of n Pomerons to m Pomerons 

g n m = n m \ n+m ~ 2 g N /2 = nm \ n+m ~ 3 g 3P /2. (7) 

No theoretical arguments or theoretical models were offered in support of this assumption, 
which certainly contradicts the pQCD approach [51 E] ■ In spite of these differences, the values 
obtained for the a tot and o~ e i are in surprisingly close agreement (sec Table 2). It is only in the 
latest version of the Durham model [3], which includes three components of the Pomeron, with 
different transverse momenta of the partons in each component (to mimic BFKL diffusion in k t ), 
that there are fairly large discrepancies in the diffractive sector i.e. a s d and add- KMR [3] find 
that at higher energies o~ s d and o~ e i have comparable values and similar energy dependence, this 
is not so in our description [2] (see Fig. 2). We note that the results presented by Poghosyan at 
this conference for a s d and Odd [IS] agree both in magnitude and energy dependence with those 
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obtained in the GLMM model. There is also disgreement in the result for the calculation of the 
Survival Probability. In [3], the Survival Probability is now multiplied by a " rcnormalizing" 
factor {(j>t)B) 2 and referred to as (S 2 ff ). The result for LHC energy is (S 2 ff ) = 0.015ioi!o5- It 
is not clear whether there is in addition a factor of (S 2 nh ) k 1/3, that needs to be incorporated. 
If affirmative, then the discrepancy between our result and that of the Durham group for (<S' 2 } 
at Tevatron energies is small, but the discrepancy becomes larger for the LHC energy range, as 
we predict that (S 2 nh ) decreases as the rapidity between the projectiles increases, while Durham 
claim little (if any) energy dependence for (S 2 nh ). 
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